function plotTemperatureField(T, r, z, time_step)
    % 获取对应时间步的温度场
    T_at_time = T(:,:,time_step);
    T_at_time1 = T_at_time / max(T_at_time(:)) * 100;  % 归一化温度场

    % 校正r和z坐标
    r_corrected = r * (2.27/0.0075) * (2.7/2);  
    z_corrected = z * (2.27/0.0075) * (2.7/2);

    % 绘制温度场
    figure;
    imagesc(r_corrected, z_corrected, T_at_time1);  % 绘制二维温度场
    axis xy;  % 设置坐标轴方向正确
    xlabel('r (mm)');
    ylabel('z (mm)');
    title(['Temperature Distribution at t = ', num2str(time_step), ' s']);
    colorbar;  % 显示颜色条

    % 移除z轴，隐藏z轴标签
    zlim([min(T_at_time(:)), max(T_at_time(:))]);  % 设置z轴的范围
    set(gca, 'ztick', []);  % 隐藏z轴刻度
    set(gca, 'zticklabel', []);  % 隐藏z轴刻度标签

    axis equal;  % 保证坐标轴比例相同
end
